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We have studied the motion of antihydrogen atoms in electric and mag- 
netic fields of arbitrary orientation. Our aim was to find an exact model 
of external (centre-of-charge) and internal dynamics which would allow us 
to describe antihydrogen motion in a realistic way. After analysing the ki- 
nematics of the charge exchange reaction producing antihydrogen and of 
antihydrogen deexcitation, we discuss the dynamics in cases of special fields 
configurations. Finally, we propose an expression for the force acting on an 
antihydrogen atom when the field configuration is arbitrary. The general for- 
mula - when reduced in the special field configurations - matches the results 
which we obtained independently before. The focus is both on the centre of 
charge dynamics and on the accurate internal quantum characteristics. 

1 Motivation and kinematics 

The motion of a charged particle in given arbitrary external electric and 
magnetic fields is governed by the Lorentz force as 

Fi = Q \Ei(R) + £ CijkXjBkiR)) , (1) 

where the i, j and k indices label the Cartesian coordinates of the vectors. 
For a system of two coupled charges (like a hydrogen atom, an antihydro- 
gen atom or a positronium), there is no such a reliable force in the given 
external fields yet, except for several special fields configurations. We have 
proposed such a force which is consistent with all crucial experimental and 
theoretical facts. Our model is simple enough to be solved analytically but 
yet sufficiently accurate (in the sense of the consistency) as well. 
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As we focus on the dynamics of antihydrogen atoms in external electric 
and magnetic fields, we suppose the fields to be arbitrary oriented at any 
point of the space and time-independent (or at least varying very slowly 
when compared to the changes in internal characteristics of an atom) . When 
we deal with the atom spectrum, we consider external fields to be completely 
constant in time (even when the atom is in motion) since the speed of inner 
atomic processes makes changes in external fields negligible; when studying 
atom as a whole, we consider all possible spatial derivatives of external 
fields. In these combined fields, we have tried to find the centre of charge 
dynamics of an individual antiatom. The particle is treated as point-like 
and its internal - quantum - properties (not necessarily independent of the 
external fields) are taken into consideration all the time. 

From the currently available technologies we consider the charge ex- 
change reaction to be the best method of the cold antihydrogen production. 
In the charge exchange, an antiproton replaces an electron in the positro- 
nium in a highly excited Rydberg state. As a product, one obtains an excited 
Rydberg antihydrogen atom: 

Ps(np s ) + p~ -> H(n) + e~ (2) 

The principal quantum numbers of a positronium and an antihydrogen are 
denoted as np s and n respectively. The kinematics of the reaction provides us 
with the initial conditions of the antihydrogen atom motion determining the 
magnitude and space orientation of the velocity vector t>g immediately after 
the charge exchange event. The kinematics of the charge exchange process 
is described by the equation [1] 

-Ug — 2 I Vpcos^g + 2^^up s cos(0p s -0g) J t>g + 

+ 2 (5* + 2 ^" f " p - ™ *^ * 0? - ds;)) = °' (3) 

where mx and vx stand for masses and velocities of the particles respectively 
and 6x = L (vx,vp)\ Rg denotes the Rydberg constant. 

By analysing (3) as a hyperplane of all possible states (vg, #g), we found 
that the angle #g does not significantly deviate from zero in any kinematics 
configuration that means the direction of its motion only slightly differs from 
original direction of the antiproton motion. 
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What remains, as far as the kinematics is concerned, is to consider the 
effect of the antihydrogen deexcitation, that is 

H(n) -)■ H(n') + 7. (4) 
We found the change in the antihydrogen velocity magnitude to be 



When an antihydrogen deexcitation event occurs, the change in the ve- 
locity magnitude of the atom is negligible. However the change of n into n' 
itself has to be taken into consideration since the dynamics in external fields 
depends on the internal state of the antihydrogen. 




2 Antihydrogen dynamics 

In our model, we assume that the dipole qualities of an antihydrogen atom 
are decisive for its macroscopic dynamics. A hydrogen atom - and there- 
fore an antihydrogen atom, too - evinces four analytically different kinds 
of dipoles: a not-induced electric dipole, a not-induced magnetic dipole, an 
induced electric dipole and an induced magnetic dipole. Moreover, one has 
to distinguish a not-induced orbital magnetic dipole and a not-induced spin 
magnetic dipole. When speaking about spin, we always refer to the spin of 
the positron; the spin of the antiproton is neglected in our model. 

We claim that a classical not-induced electric dipole d (for instance a 
water molecule) placed in a stationary external electric field E = E(R) 
experiences a force expressed in the Cartesian coordinates as 

where d stands for the projection of the dipole vector d into the direction 
given by the local vector E. When the quantum nature of an antihydrogen 
atom is to be effectively implemented, the d quantity is the function of the 
proper quantum numbers. 

The force acting on a dipole can be derived using two approaches. The 
first method treats the dipole as an object of a finite size - this corre- 
sponds to the separation of two charges defining the dipole. The system has 
more than three degrees of freedom which may be however reduced: The 
torque equation N = dx E prescribes that E is the local preferred direction 



3 



around which the d vector undergoes internal motion with a well-defined 
time-averaged value d. This holds for time-independent external fields. The 
force on this particularly defined dipole is therefore fully responsible for the 
dipole's macroscopic motion and equals (6). This approach was inspired by 
outcomes in [2] and by unpublished result [3]. 

The second possibility is to take the classical expression W = —dE as 
the potential energy of a not-induced electric dipole in an external electric 
field (with the quantum conditions still easily applicable on d) [4]. In this 
case, the force is 
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The situation for a not-induced magnetic dipole p placed in an external 
magnetic field B = B(R) is similar: Its potential energy W = —jlB leads to 
the same force on the dipole 



_ A Bj dBj 



(8) 



as it is obtained by solving the force on an elementary current loop. The 
p, stands for the projection of \x into that preferred B direction. Moreover, 
the internal dynamics is analytically solvable and known as the Larmor 
precession. 

For induced electric and magnetic dipoles separately, the dynamics is 
very similar: The forces resemble (6) and (8), only the dipoles are explicitly 
proportional to the field magnitude with accordance to their definition. In 
these cases we use olE instead of d in (6) and (3B instead of p, in (8) re- 
spectively. 

The result (6) is also applicable to the antihydrogen atom in such a field 
if d is substituted with a proper function of the internal quantum numbers 
of the antihydrogen. The shifts of spectral lines of the antihydrogen atom in 
an electric field [5] equal 



W = -a enpE, 



(9) 
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where ao denotes the Bohr radius, e is the elementary charge, n is the prin- 
cipal quantum number of the atom and p stands for the parabolic quantum 
number (an element of the {— n + 1, —n + 2, n — 2,n — 1} set). If compared 
to W = —dE, one has d = — ^a^enp and the force on such a particular anti- 
atom in an external electric field is 

_ 3 En QEa , . 

F ' = -2 W £iM? (10) 

which is (6) in fact, with the quantum character of the atom implemented. 

The other three dipole cases for the particular antiatom can be evaluated 
in very similar way. Namely for that of a not-induced magnetic dipole in a 
magnetic field, the energy (Zeeman effect) shifts [6] are 

mB, (11) 

2m e ' v ' 

where m e is the rest mass of an electron, h is the reduced Planck constant and 
m stands for the magnetic quantum number (for a given n it is an element 
of the same set as the parabolic quantum number is). When compared to 
W = —ptB, one obtains 

eh ^BjdBj 



2m e B dX t 



when [i = —^-m is substituted into (8). Note that both (11) and (12) hold 
only if the spin of the electron (positron) is neglected - the effect of spin 
will be discussed later. 

The antihydrogen spectra themselves prove to be reliable starting point 
for treating the atom dynamics. The question now is, how these energy shifts 
look like when both electric and magnetic fields are present (that is which 
expression reduces to the (9) and (11) energies in the single-field limits). The 
expression for W should be explicitly dependent on the angle 7 between the 
two local field orientations (7 = l(E,B)) [7]. 

The most straightforward way to construct the antihydrogen spectra 
would be a simple sum of (9) and (11) - the Stark and Zeeman energies. 
However, such an object, described by 

W = --a enpE - —mB, (13) 
2 2m P 
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does not satisfy the condition of the angular dependence and is therefore 
incorrect. The accurate expression has to be found - this is the core of this 
paper. 

With the spin still not taken into account, for any local electromagne- 
tic field, such a Cartesian coordinate system can be found that the most 
basic Hamiltonian [8] of a hydrogen atom in the simultaneous electric and 
magnetic fields reads 

- 2 2 fc 

P 6 Gil 

H = H BL3 — eE (x sin 7 + z cos 7) . (14) 

2m e Aireo? 2m e 

The perturbation theory leads to (9) when B = is set in (14) and we get 
(11) for E = as well. Unfortunately, there is no known way for a general 
fields configuration how to analytically evaluate the spectrum of the matrix 
whose elements are given by the integral 



w (i'm')(im) = J R * d/Y * m , (+JL- B L 3 - eEx^ R n iY lm dV, (15) 

wh ere (j£ - ^) (R nl Y lm ) = (R nl Y lm ). 

There are alternative methods of finding the spectra, namely the Som- 
merfeld quantum conditions which give exactly the same results as (9) and 
(11). What makes them important is that the combined problem has been 
already solved using them, as described by Max Born in [9], for instance. 
The final spectrum reads 
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(16) 



where two new quantum numbers tia and ub (instead of p and m) were 
introduced. Both ua and tib are independently elements of the set 

f n—l n—1 _i_ 1 1 n—1 i i n— 1 \ 

\ 2 ' 2 ' ' " ' 2 ' 2 J" 

Due to the fact that (16) does not reflect the effect of the spin of the 
electron, the outcome has to be slightly modified. It seems impossible to 
implement the effect generally, even in the absence of the electric field. 
There are in fact three distinct not-induced dipoles in a hydrogen atom: not- 
induced electric dipole, not-induced orbital magnetic dipole and not-induced 
spin magnetic moment which have to be distinguished. The expression (16) 
connects the first two together with no limitation on the field strengths. 



On the contrary, the connection of the two not-induced magnetic dipoles 
depends on the field strength [8]. However, to maintain the ambition of a 
universal equation of motion, we suggest the spin effect in a way we consi- 
der simple and minimally harming the accuracy - by adding the term f^sB 
(where s = ±5) to (16). Moreover, when the absolute values are expanded, 
also the explicit angular dependence on 7 becomes 

1/3 \ 2 f eh \ 2 3a e 2 nh eh , . 

+ ^- W nE) +(— B) -^—EB CaS y„ B + -sB. (17) 

We found that 

-^—BL, eEr 

This relation also implies that the associated quantum numbers p and m are 
well-defined only when commutator (18) equals zero, that is in the pure Stark 
and Zeeman cases, for instance [10]. Therefore the new quantum numbers 
riA and ns were introduced instead of p and m (in our model, the spin 
quantum number s is well-defined for any fields configuration). The total 
spin-less degeneracy of the ra-th manifold equals 

n— 1 +1 n— 1 n—1 n— 1 / -i s 

E E i = E(^ + 1 ) = 2 E^ + E 1 = 2 - !ii! ^ 12 + - = ^ (is) 

1=0 m=-l 1=0 1=0 1=0 

which has to be satisfied for the new quantum numbers as well and 

I n—1 I n—1 

2 "r 2 

E E i 

„ n—1 „ n—1 

n A - — n B - — 

The Stark and Zeeman forces are two special cases of the (16) or (17) 
results. Note that instead of (12), where the effect of the spin was not taken 
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into account, we obtain 



„ eh , . Bj dBj 

F -=-s: ( "' +2 ^iw (21) 

Aside of the Stark and Zeeman effects, we found two more special cases. The 
first is the case when the electric and magnetic fields are parallel - according 
to (18), the two perturbations are independent and the "old" quantum num- 
bers p and m are well-defined. The potential energy in this case is therefore 
the simple sum (13) modified by the term ^sB responsible for the spin 
effect. Hence the force on an atom in parallel fields is the negative gradient 
of that potential energy and reads as 

eh , . Bj dBj 3 Ej dEj . . 

r ' = -H (m+2 «iiw-^" p ^i«' <22) 

6 J=l % 3=1 1 

The final case which we were able to solve independently of (17) is such 
configuration of fields when the magnetic field is much stronger than the 
electric field. We analysed this case using the fact that the magnetic field 
is dominant and therefore we constructed the potential energy as W = 
2^- (m + 2s) B + ^aoenpE cos 7 - the angle dependence is present at the 

electric term since the orientation of the B vector is significant for the atom 
and the second term as whole is considered to be a small correction. The 
force on an atom in such fields is 



eh , . -B,- dBj 3 JU Ej dEj , . 

"■ = -H (m+2s) ^i^- 2 w ^ ew, m ->- (23) 

e j=i 1 j=i 1 

Now let us return to the general case for which the potential energy 
is given by (17). In very strong fields, the effects of the induced dipoles 
(quadratic effects) have to be taken into account. We assume that the two 
terms —^aE 2 and — \fiB 2 responsible for them can be simply added to 
the potential energy since their effects are independent of the not-induced 
dipoles in the first approximation. The polarisabilities a and f3 depend in 
a bit more complicated way on ha, ns and n. The total external force (as 
the negative gradient of the potential energy) on such a model of an atom 
therefore is: 
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m e ^ B dXi ^ 3 dXi H ^ 3 dXi v ' 

What makes us confident about this result is that all the independently 
solvable cases - the outcomes (10), (21), (22) and (23) - are exactly the limits 
of (24) for the particular fields configurations. In the four independently 
solvable cases, we found that (and how) the quantum numbers p and m can 
be linked to ua and rig. 

All the terms are proportional to the field derivatives hence there is 
no force on an atom in homogeneous fields. All the weak field terms are 
antisymmetric in ua, ub and s. This allows using the formula (24) for both 
matter and antimatter with only the redefinition of these quantum numbers. 



3 Conclusion 

We have analysed the spectrum of an antihydrogen atom in time- independent 
external electric and magnetic fields of arbitrary mutual orientation. We have 
discussed the kinematics of a possible deexcitation and the kinematics of a 
charge exchange reaction - which we consider the best method for the cold 
antihydrogen production. Charge exchange kinematics determines the initial 
conditions for the motion of the antihydrogen. 

In the search for the antihydrogen dynamics itself, we studied the ele- 
mentary and analytically solvable cases, such as an electric dipole in an 
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external electric field for which the Stark-like hydrogen is a special case. In 
this sense we have justified the spectral approach to be reliable and equi- 
valent to the potential energy since the antihydrogen possesses three degrees 
of freedom in our model. We used the result derived by Max Born [9] to be 
such an accurate spectral term for an antihydrogen atom and we implemen- 
ted the effect of a spin and possible strong-field effects. The potential energy 
satisfies all the experimental properties and so does the external force acting 
on the antiatom, which is (24). 

We hope that the equation will prove itself useful and it will help to 
design experimental devices and carry out measurements on antihydrogen 
atoms. We are especially interested in the gravitational properties of an- 
timatter and we hope that the theoretical understanding of the antiatom 
motion in the combined fields will contribute to the precise understanding 
of systematic errors in measurements of the atihydrogen in free fall after its 
acceleration. If we succeed in effective trapping of antiatoms, much more 
experiments would be possible. The same holds for positronium. 
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